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Abstract 

We study a many-server queuing system with general service time distribution and state 
dependent service rates. The dynamics of the system are modeled using measure valued 
processes which keep track of the residual service times. Under suitable conditions, we 
prove the existence of a unique fluid limit. 
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1 Introduction 

In recent days, many-server queuing systems have received much attention due to its appli- 
cations to call center. Thus it has become important to study its asymptotic properties to 
gain insight into the behavior of these systems. Studying different scaling limits (fluid or diffu- 
sion scaling) are established tradition in queuing theory. In the celebrated work of Halfin and 
Whitt [TJ, it was shown that for Poissonian arrival and exponential service time, with positive 
probability, there is a positive queue in the asymptotic regime if the arrival rate A n and the 
number of servers n both goes to infinity in a manner that X n = n — (3y/n for some /3 > 0. 
Fluid and diffusion limits for the total number of customers in a network with time varying 
Poissonian arrival and stuffing was obtained in |12| . This work was generalized in for 
Gt/Mt/st queuing networks with abandonment where the authors studied long-time behavior 
of the fluid limits. 



A recent statistical study by Brown et. al. [2] suggests that it may be more appropriate to 
consider non-exponential service times. In particular, it is log-normal in some cases as shown 
in [2]. This emphasizes the need to consider many-server model with generally distributed 
service times. In |15j . Whitt considered G/G/n network with abandonment and proposed a 
deterministic fluid approximation. The author proved convergence for discrete time model. In 
[TU] , Kaspi and Ramanan considered G/G/n model and obtained a measure-space valued fluid 
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limit. Later Kang and Ramanan generalized this work by allowing the customers abandonment 
in [8]. In [I], Atar et al. studied multi-class many-server queues with fixed priority and 
established the existence of unique fluid limits. Kang and Ramanan studied ergodic properties 
of the scaled processes for GI/G/n+G model and its relation with invariant states of the fluid 
limit in [9]. Reed in [T3], established the fluid and diffusion limits of the customer-count 
processes for many-server queuing system under the finite first moment assumption on service 
time distribution. In [IT], Zhang obtained the fluid limits for GI/G/n+G queuing systems. 

All of the above models consider servers that serve the customers at a constant rate 1. 
In this work, we allow the servers to adjust their service rate depending on the number of 
customers in the system (equivalently, the number of customers in the queue). It is often 
useful to increase the service rates when the queue length is large. Management may also be 
interested to adjust the service rate depending on customers feedback. State dependent arrival 
and service rate were first introduced in [16] for conventional heavy traffic approximations. 
In case of single server models, processor sharing model is an example where service rate at 
any instant of time depends on the number of customers in the system. For some recent 
developments on processor sharing, we refer [6], [18], [13]. 

In this work, we consider a system with n-homogeneous servers. Customers arrival is 
given by a renewal process and customers are served under FCFS policy. Arrived customers 
do not leave the system until served. Let X™ denote the number of customers in the n-th 
system at time t > 0. Define X" = ^X™. The service rate of each server at time t is given 
by k n (Xt) for some bounded map k n on [0, oo). Note that, k n {X^) could be for non- 
zero X™. The system is described by (Q2,Z™) where denotes the number of customers 
waiting in the queue at time t and Z™ is a non-negative Borel measure on (0, oo) such that 
Z^iC) denotes the number of customers in service with their remaining service requirements 
in C, for C 6 £>((0, oo)). Thus Z n is a measure- valued process that keeps track of customers 
remaining service requirements. Measure- valued processes that keep track of residual service 
requirements of individual customers, have been considered earlier in literature (see [3], [6], 
|18j . |17j). In [IH], the authors used measure- valued processes that keep track of the time spent 
by the individual customers in service. Also their proof relies on the fact that there exists a 
compensator for the departure processes (see Corollary 5.5 there). Since we are allowing our 
service rate to be dynamic depending on X n , getting an explicit compensator for analogous 
processes as in jlO] is a hard problem. In this work, we closely follow the approach in [18] (see 
also [IT])- We show that the fluid limit of -rZ^ is uniquely determined by an integral equation, 
referred as fluid model equation (([9]) below). A similar type of equation was also obtained in 

m- 

This paper is organized as follows. In Section 2, we introduce our GI/G/n model with our 
basic assumptions and state our main result. The uniqueness part of the main result is done 
in Section 12.11 In Section 3, we prove relative compactness of the stochastic pre-limits and 
characterizations of the limits are done in Section 4. Finally, in Appendix we prove existence 
result for solution to certain integral equation and recall some results from [18] those are used 
in this paper. 
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1.1 Notations 



The following notations will be used throughout this paper. By N,R, we denote the set of 
natural numbers and the set of real numbers, respectively. Given a, b G M, the maximum 
(minimum) is denoted by a V b(a A b). We use a + for a V 0. We define = [0, oo). For any 
A C [0, oo), we define A € = {x > : inf ag ^ \ x — a\ < e}. For any the sets (x, oo), [x, oo) 

will be denoted by C X ,C X respectively. For any topological space S, Cb(S) denotes the set 
of all real valued bounded, continuous map on S and B(S) is used to denote the Borel a- 
field of S. C([a,b],S) will denote the set of all continuous function from [a,b] to S. For any 
/ G C{[a,b],R), we define \f\ ab = sup xgM \f{x)\. 

The set of all non- negative finite Borel measures on [0,oo) is denoted by M and 
denotes the subset of Ai containing all the measures having no atom at {0}. For any fj, G Ai 
and Borel measurable function g on [0,oo), we define {g,fj) = J gdfj,. For /ii,/i2 G -M, the 
Prohorov metric is defined by 

p(/ii,/i 2 ) = inf{e > : m(A) < /J. 2 (^4 e ) + e, fJ, 2 (A) < m(A 6 ) + e,for all closed A C [0,oo)}. 

It is well known that (A4,p) is a Polish space (see Appendix in [3]). Also this topology is 
equivalent to the weak topology on Ai which is characterized as follows: /U n — > /i in weak 
topology if and only if 

(/,/in)^(/,M>for all/GC 6 ([0,oo). 

Given any Polish space (E,ir), D([0,oo), E) denote the space of functions that are right- 
continuous with finite left limits (RCLL). Endow the space D([0, oo), E) with the Skorohod- 
Prohorov-Lindvall metric or J\ metric, defined as 

d(&0')= pf (||/||°V I™ e- u d u (<i),<i)',f)du), 0,0' eD([0,oo),E) 
f eT v Jo J 

where 

d u ((f>, (/)', f) = sup[vr(0(i A u), 4>'(f(t) A u)) A 1], 
t>0 

and T is the set of strictly increasing, Lipschitz continuous functions from [0, oo) onto itself, 
with 

f(t) - m 



sup 

0<s<t 



log' 



t - s 



< oo. 



As is well known [5], D([0,oo), E) is a Polish space under the induced topology. 
We use " =^ " to denote the convergence in the sense of distribution. 



2 Queuing model 

In this section, we describe our GI/G/n model and the measure valued state descriptors. We 
assume that for each n, all the stochastic variables below, are defined on probability space 
(f2 n , J- n , P n ). The system contains n identical servers. Each arriving customer has a single 
service requirement and is served by a single server. The customers are served by FCFS policy 
and they leave the system once their service is completed. We do not allow the customers to 
renege the system until their job is done. We also assume that the system works under work 
conserving policy i.e., all the servers are busy if there is a queue. We assume the following: 
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• Customers arrive according to a renewal process E™ with mean inter-arrival time for 
some A™ > 0. 

• Service requirement of the arriving i-ih. customer is given by vf where {w"}^^f is an 
positive valued i.i.d. sequence with common distribution u 11 . 

Upon arrival the customers join the queue if all the servers are busy. At time t, all the 
servers serve the customers at a rate k n {Xf) where Xf = — - and Xf denotes the number of 
customers in system at time t. It is reasonable to assume that k n : M+ — > M+ is a bounded, 
Borel measurable function. Let r", i S N, be the time when the i-th customer starts its service. 
Then for t > 0, the remaining service of i-th customer is given vf — J „ k n (X™)ds (non positive 

4 

quantity implies that customer's job is completed) provided r™ < t. We use negative indices 
to denote the customers in system at time t = 0. Let Xq denote the number customers at 
time t = with remaining job vf, i = —Xq + 1, . . . , — [Xq — n] + , for the customers in service 
at time t = where vf,i = —Xq + 1, . . . , — [Xq — are some random variables defined on 
(f2 n , Tni IPn)- Also let Qq denote the number of customers in the queue at time t = 0. Hence 
Q n Q = [XZ-n]+. 

For t > 0, denotes a measure in Ai + such that Z^{C) denotes the number of customers 
in service with remaining service requirement in C for C C ((0, oo)). Hence the total number 
of customers in service at time t is given by Z™ = Z™(0, oo). Let be the number of 
customers waiting in the queue. Define = — Qf. Then + 1 denotes the index of 
the head of the customers in the queue waiting to be served. For < s < t < T, we define 
S n (s,t) = f s k n (X™)du. Hence a precise description of Z n is given by 

-Qo B t 

znc)= Yl tv?(c+s n (o,t))+ Yl ^(c+s n (r?,t)), a) 

i=-X$+l i=-Q" + l 

for C S i3((0, oo)). Additional obvious relation satisfied by the processes, for t > 0, are as 
follows: 

AT = Q n t + Z?, (2) 

Q n t = [X? - n}+. (3) 

It is easy to see from (jj)) and ((3j) that = A" A n. We extend Z™ on B([0, oo)) by setting 
Z™({0}) = 0. It is easy to see that Z n takes values in D([0, oo), M). 

2.1 The fluid model 

In this section, we define the fluid model and state our main theorem. We also state the set of 
assumptions that are used to prove this result. We define 

yn i rya yn 

y-n _ *H _ 1 jrn f\n _ yn _ 

n n n n 
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The fluid scaling of arrival process Eg is define as Eg = —Eg. Thus the fluid pre-limit equations 
are (in analogy with dl])-©) given by 



Zt{C) = - jr Se? (C + S n (0,t)) + - jr 5 V? (C + S n (rr,t)), (4) 
for C € B([0,oo)). 



n * — ' 1 n 

«=-Qo+ 1 



Xg = Q? + Zg, (5) 
Ql = [Xg - 1]+. (6) 
We assume the following conditions: 

Condition 2.1 (a) — -> A /or some A E [0, oo) and =>• At in t/te sense of distribution 
in D([0,oo),R + ); n 

(b) There exists a probability measure v with bounded, Lipschitz continuous, density g : R+ — > 
K+ such that v n — > v as n — >■ oo; 

(c) (Zq,Qq) =>■ (^OjQo) in A4 x M + as n — >■ oo where the function F(x) = Zq([x, oo)) is 
Lipschitz continuous on and (-Zq,Qo) * s a deterministic element in Ai x R + . 

Since At is a continuous, deterministic path, one obtains the convergence of the scaled arrival 
process in probability i.e., for any T, e > 

lim P n ( sup \Eg - \t\ > e) = 0. (7) 

rwoo o<t<T 

We impose the following condition on the state dependent service rates: 

Condition 2.2 There exists a bounded, Lipschitz continuous map k : R+ — > M + such that 
k n —> k as n —> oo uniformly on compact subsets o/R+. 

Let G(-) be the distribution function of v i.e., G(x) = u([0,x]). Define G c (x) = 1 — G(x). Also 
from ([5]) , © and Condition 12.1( c) the followings hold: 

(1) ^o) = %o, Xo = Qo + Zo, Qo = [-Xb - 

Theorem 2.1 Assume Conditions \2.1l POl to /io/d. T/ien as n — > oo, (Z n ,Q n ) =>■ (Z,Q) in 
D([0, oo),Ai x R + ) where (Z,Q) is uniquely defined by the followings: 

• /or a// t > 0, 

Z t ((0, oo)) = Z t , X t = Q t + Z t , Q t = [X t - 1] + , (8) 

• for all t > and x > 0, 

£*(£*) = F(x + 5(0, t)) + / G c (x + 5(«, t))dB s , (9) 

J o 

where S(s, t) = J s * k(X u )du and B t = At — Qt- We re/er ([9]) as /Zuid model equation and 
will be denoted by (k, A, z/). 
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Proof: The existence of a limit satisfying the above properties will be proved in Section 3 and 
4. So it is enough to prove the uniqueness of the limit here. Since B™ is nondecreasing and 
we can have pointwise convergence for each subsequential limit (by Skorohod representation 
theorem), Bt is also non-decreasing. Therefore ([9]) makes sense. Also for any a > 0, 

Z t ([0,a))<O(a), 

implying Z t ({0}) = for all t > 0. Now putting x = in ([9]), we hav^j] 

Z t (C ) = F(S(0,t)) + [ G c (S(s,t))dB s 

Jo 

= F(S(0,t))+x[ G c (S(s,t))ds - [ G c (S(s,t))dQ s 
Jo Jo 

= F(S(0,t))+x[ G c (S(s,t))ds-Q t G c (0) + Q G c (S(0,t))- [ g(S(s,t))Q s d(S(s,t)), 
Jo Jo 

where we used integration- by-parts formula in the last line. Since G c (0) = 1, using ([8]) we have 

X t = Z t + Q t = F(S(0,t)) + Q G c (S(0,t)) + \ f G c (S{s,t))ds 

Jo 

" (X s -l)+g(S(s,t))d(S(s,t)). 



By our assumptions on F(-),G(-) and k(-), we see that S(-,t) is Lipschitz continuous on [0, t] 
and hence Xt is continuous in t. Therefore 

X t = Z t + Qt = F(S(0,t))+Q G c (S(0,t)) + X I G c (S(s,t))ds 

Jo 

+ [ (X s -l) + k(X s )g(S(s,t))ds. 
Jo 

Hence by (jH]), Z t ,Qt are continuous in t. Now defining H\{x) = F{x) + QqG c {x),H2(x) = 
XG c (x), H±{x) = g(x), H$(x) = k(x) for x > and extending these maps on (— oo,0] by their 
respective values at 0, we have 



X t = H 1 {S{0,t)) + f H 2 (S(s,t))ds+ [\x s -l)+H 4 (X s )H 5 {S{s,t))ds. 
Jo Jo 



(10) 



By Lemma 15.91 in Appendix A, Xt is uniquely defined on [0, T] for all T > 0. Since Qt = 
[X t — 1] + , Qt and hence B t is unique. Therefore ([9]) implies that Z t is uniquely defined on 
C x . Since {C x ,x £ M + } defines uniquely a Borel measure on (R+, B([0, oo))), Z t is uniquely 
defined by ©. Hence (Z, Q) is unique in D([0, oo), M x R + ). □ 

Remark 2.1 One can relax the conditions on G(-) depending on the properties of service rate 
k(-). For example, if k(-) is constant then it is enough to impose continuity on G(-) instead of 
Condition\EM h ) J731 \nj). 



corrected with - below, last line 
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3 Tightness of the pre-limit processes 

In this section, we study the compactness properties of the pre-limit processes. From (j3|), we 
get the following equation 



Zl\C) = Z™{C + S n {s,t)) + ±- J2 6 v n(C + S n (r?,t)), (11) 

i=B»+l 

for all < s < t and C G S([0,oo)). Define E n (s,t) = E? - E%. From ©, it is easy to see 
that given T, e > 0, 

lim P n ( sup \E n {s,t) - X(t-s)\ < e) = 1. (12) 

n->oo 0<s<t<T 



Defining B™ = — , we have 



B? = E?-Ql (13) 



We recall the following characterization of compact subsets of M in Prohorov topology from 
[3] (Theorem A2.4.I). 

Definition 3.1 A set K C M is relatively compact if and only if sup„ eK < oo and 

there exists a sequence of nested compact sets Cj C M+ such that L)Cj = M + and 

lim sup ^(C 7 C ) = 0. 
The proof of the following lemma is same as [181 Lemma 5.1]. 

Lemma 3.1 Fix T > 0. There exists a sequence {e£;(n)} such that e^(n) — > asn->oo and 



( sup \E n (s,t) - X(t- s)\ <e E {n)) > I - e E {n) 

0<s<t<T 



for all 77, G N. 



We define £l E = {sup 0<s<t<T \E n (s,t) — X(t — s)\ < e^(n)}. The following lemma proves 
compact containment properties of (2 n ,Q n ). 

Lemma 3.2 Assume Condition \2.1\ to hold. Fix T > 0. Then for any positive rj there exists 
a compact set K C M and K > such that 

liminf F n (Z? G K and Q? <K for all t G [0, Xl) > 1 — r/. 

n— >oo 

Proof: By Condition 12.1( c). there exists a positive integer Mq such that 

su P P n (Q£ > M ) < ^ 

n o 

Since Q? < Q% + E%, choosing K = M + 2 AT we get from (JT5J) that 

limsupP n (Q? > K for all i G [0,T]) < 5 (14) 

n— >oo 4 
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Define i?£ = ± E i= Lgn +1 o"„« G X. From (HU), we note that 

Z?(.C x )<ZS(C x ) + 0?(C x ), (15) 

for all x G R+ and i > 0. For m G Z,i > 0, define £ n (m,£)_= £ ES+S Reca11 
the definition of f^(M, L) (see (|B2|) in Appendix) and function / from Appendix B. Define 
^i = {Qo — M)}- Therefore using Lemma 13. II we have for all large n, 

F n (ni nn?n fi£(M + 1, *Q) > 1 - 1. 

Choose n large enough so that e^(n) < 1. Then on fi^ n n f2^(Mo + 1, K) we have for all 
large n 

(/, 0») < z/ 1 ) + 1 < (/, P) + 1 < Mi, (16) 
for some positive constant Mi. Hence using Markov's inequality and (fT6|) we get 

0?(C X ) < j^-M 1 , (17) 

on S]|nOJn ^^(M + 1, A") for all t G [0, T] and all large. Again by Condition EjTc), there 
exists a compact Ko C M such that for all large n 

F n (ZZ G Ko) > 1 - |. (18) 

We denote the above event by O^- By Definition 13 .1\ there exists a bounded function g : R+ — >• 
such that lrm a; _ > . clo g(x) = and on Q^i 

2q{R+) < p{0), Z%{C X ) < q{x) for all x G R+. (19) 

We define 

K = {/i G : < p(0) + K, fi(C x ) < g(x) + -=— -Mi V x G 

/(«) 

By the property of g and /, K is a compact subset of Ai. From (|15p . (|17p and (|19|) we see 
that for all large n, on 0| n fif n O^(M + 1, if ) D 

g k, 

for all t G [0, T]. Also for all large n, P n (0£ n fi? n nj(M + 1, iT) D fig) > 1 - Hence the 
proof follows from (|14p . □ 



Lemma 3.3 Assume Conditions \2.1l \2.2\ to hold. Fix T > 0. T/ien /or each e,n > 0, f/iere 
exists a k > suc/i i/mf 

liminf P n (sup sup Z t n ([x, x + «]) < e) > 1 - r/. (20) 
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Proof: Using Condition 12.1( c). one can prove that for any e, n > there exists a positive k 
such that 

liminfP n ( sup ZZ([x,x + k}) < e/2) > 1- -. (21) 

In fact, the proof is same as the proof of (79) in [16| . We denote the above event by and 
the event in Lemma by J) J. Define ftg = 0£ n n ^ n fi£(M, L) for L = 2AT + K and 
M = [L\ + 1 . From Lemma 13.21 it is easy to see that there exists K > such that 

liminfP n (n?) > 1 — 77. 

n— >oo 

From pip , we see that for any k > and t G [0, T] 

Z t n ([x,x + K ]) = J n ([x,x + K]+5 n (0,t)) + - V ^n([x,x + K ] + 5"(rf,t)). (22) 

n z — ' 1 
«=-Qo+i 

On Q3, we have sup x6]K+ ([x, x + k]) < e/2. Hence choosing x = x(uS) = x + S n (0,t) we have 
on Jig , supj-gjj Zq([x, x + k] + ,S n (0, t)) < e/2. Thus we need to estimate the second term on 
O". So we denote the second term by E t . 

For any 5 > 0, we consider a partition = to < t\ < ■ ■ ■ < t r = t of [0, t] with \U + x — ij| < 5 
for k = 0, 1, 2, . . . , r — 1. Since — = #0 > we nave 

1 -Bf" 

7" — 1 fe+1 

H ' = E^ E 5v?([x,x + K]+S n (T?,t)). 

k=0 i=B? +1 

From (jB2|) . on fi^(M, L), we have for all fc = 0, 1, 2, . . . , r - 1, 

-Qo < ^ < E?, and so, < B^ - B% < E? + QJ, 

and 

max sup sup|(/,£ n (m,£)) | < e A (n). 

Hence for all m G (— nM, nM),£ G [0, L] and for all a, 6 G M + , a < 6, we have 

(X [o ,&],£(fM)) <% M ,^) + 2e A (n), (23) 

on^(M,L). Nowfortfc < rf < t k+1 , [x, x + k] + 5 n (r™, t) C [x + S n (i fc+1 , t), x + K + S n (£ fe , t)]. 
Now fixing a = x + S n (t k+ i,t),b = x + k + S n (t k ,t) and observing that, on $7^ n f^, S ti G 
(—nM, nM) (for above choice of M) and — ^ G [0, L] for all n large , we have 

n+ i 

- Y, fy([x,x+K]+S n (T?,t)) < (B? i+1 -BZ)v n ([x + S n (t k+1 ,t),x + K+S n (t k ,t)]) + 2e A (n), 

for k = 0, 1, 2, . . . , 1 — 1. Since — > v in Prohorov metric, for any ei > 0, there exists no G N 
such that for all n > no, and closed set C C R+ (see Notations) 

v r \C) <v{C e ^) + e x . 
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Hence combining above two we have 

Bf 

z i-\-l 

- ]T S K ([x,x + K] + S n (rr,t)) < (B? i+1 -BZ)v([x + S n (t k+1 ,t)-e l ,x + K + S n (t k ,t) + 

i=BV-+l 

+(Bl +i -Bl)e 1 + 2e A (n), 

on S72 n for all n large. At this point, we observe that sup se [ ,r] Xs ^ su Psg[o,t] + 1 < 
L+l on for all large n. Hence by Condition ^. 2| on sup[ T ] fc n (X") < sup[ L+1 ] /c n (x) < 
su P[o,l+i] k(x) + 1 < M2 for some positive constant M2 and for all n large. Hence \S n (t k ,t) — 
S n (tk+±, t)\ < 5M2 on for all n large. By Condition I2.1f b). we can choose 5 and k small 
enough so that on 

u([x + S n (t k+1 ,t) - ei,x + k + S n (t k ,t) + ei]) < 

for all ei small enough and all n large. Hence summing up the above expression we have for 
all t £ [0, T] 

3* < L.^- + Le 1 + 2re A {n), 

on fig for all n large. Since ei, €A(n) do not depend on r and x, we can choose them small to 
make the right hand side smaller than | for all n large and x £ IR + ,i £ [0,T]. The proof is 
done from (|22|) and definition of fi??. □ 

For any path (f> £ -D([0, 00), E) where (E,ir) is polish space, the 5-oscillation of 4> on 
[0,T], T > 0, is defined as follows 

w(4>,5,T)= sup 7T(0(s),0(t)). 

s,te[0,T],|s-i|<<5 

The following lemma gives the oscillation bounds on the stochastic process 

Lemma 3.4 Assume Conditions \2.1l \2.S\ to hold. Fix T > 0. Then for each e, r] > 0, there 
exists 5 > such that 

liminf F n (w(Z n ,S,T) < e) > 1 - 77, 

n— >oo 

liminfP n (-u;(Q n ,<5,T) < e) > 1 - 77. 

n— >oo 

Proof: Let t,s £ [0, T],s < t and |i — s\ < 5. Let -D™ be the number of customers finished 
their job by time t. Then it is easy to see that 

Bf 

D?-D n s <Zl l %S n (s,t)]+ ]T 5 v? ([0,S n (s,t)])- 

i=B"+l 

Recall the event Q§ from Lemma 13.31 By definition, on fi^, we have -B™ £ (—nM,nM) and 
B™ — B™ £ [0,L] for all n large. Since sup[ T ] X r s l < sup[ T ] Q r s l + 1 < L + 1, we can choose 
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5 small enough so that \S n {s,t)\ < K\ on fig for all n large where K\ = K\{§) — > as S — > 0. 
Therefore by the definition of fi^(M, L) and ()23[) we have 

i ^ S vf ([0,S n (s,t)})<L( X[o , Kl] ,v n )+*A(n), 
i=By-+i 

on fig for all n large. Since v n — > v, by the same reasoning as in Lemma 13.31 we can choose 8 
small so that 

1 £ 5 v n([0,S n (s,t)])<e, 

on fig for all n large. If we denote the event in (|20p by fig , then on fig n fig 

Z^[0,S n (s,t)]< 4«[0 )Kl ]<e 

for all n large and 5 chosen small enough. Hence with this choice of 5, P n (fig n fig) > 1 — 2r] 
and 

i(A ra -^)<26 
n 

for all n large. Since = X £ + £ t n - , we have |X s n - | < |£f - + I (D t n - L>") < 3e 
on fig n fi^ for all n large provided 5 chosen small enough. Therefore with this choice of 5, we 
have (using ([6])) 

io:-on<3e 

on fig n fig for all n large. Hence the second claim follows by replacing e, 77 with e/3, r//2 
respectively. 

Now we prove the first claim. We note that for any t,s€ [0, T], \B% - B?\ < \E% - Ef\ + 
\Q n s -Qt \ (from $13])). Hence on fi^nfig 1 , |-B"-5 t n | < 4e for all n large provided |i-a| < 6A^. 
Let C C M + be closed and C 6 be its e-enlargement. Now choose 5 small enough so that n\ < e 
and so C + S n (s, t) C C Kl C C e on fig for all n large. Hence from (jlip we have 

(C) - J"(C e ) < J"(C7 + S n (s, t)) - Z n s {C^) + \B n s - B»| < 4e, 

on fig n fig for all n large. Again for any c G C and <S n (s, i) < «i, we have dist (c—S n (s, t),C) < 
e implying c G C e + S n (s, t) and so C C C e + 5 n (s, t). Hence from (JTTD 

Z n s (C) - Z?(C e ) < Z:{C) - Z2(C + S n (s,t)) < 0, 

on fig n fig for all n large. Hence for all closed set C G B([0, 00)) we have 

Z n s (C) < Z?{C ie ) + 4e and Z?{C) < Z^{C ie ) + 4e, 

on fig n fig for all n large. Hence p(Z™,Z™) < 4e on fig n fig for all n large. Thus first claim 
follows by replacing e, r] with e/4, 77/2 respectively. □ 

Now we introduce a weaker oscillation function w' on D([0,oo),M x Define metric 

d'(-,-) = max{p(-,-),| • |} on Ai x M + which induces a separable complete metric on it. For 
any t/j G D([0,oo),M x M+) and T, 5 > define 

w/^, 5, T) = inf max sup d! (ip(s),ip(t)), 
li 1 s,te[ti-i,ti) 
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where {tj} ranges over all partition of the form = to < t\ < . . . < tj = T with mini<j<j(tj — 
ti-l) > 5 and j > 1. It is easy to see that for 5 > we can have a partition ij of [0, T] such that 
mini<j<j(tj — ti-x) > 5 and maxi<j<j(ij — ii-i) < 2<5 and hence for any ^ G D([0, oo), xIR_|_) 
we have 

w'(tp,5,T) < w(tf;,28,T). (24) 
Hence from Lemma 13.41 and (|24p . we get that for any T,e,rj > 0, there exists 5 > such that 

limmiFJ w'((Z n ,Q n ),5,T) < e) > 1 - rj. (25) 

n— >oo 

For ip G L»([0,oo),X x E+), define 

POD 

J(V>) = / e _s [J(^, s) A l]ds, where J(<p,t) = sup d'(ip(s), ip(s-)). 

JO 0<s<t 

Again it is easy to see that J(tp, t) < w(ij), 5, t) for any 5 > and hence J(ip) < u>(^>, 5, T) + e~ T 
for all T, 5 > 0. Thus applying Lemma |3.4| we get that for any e, rj > 0, 

liminf P n (J((Z n ,Q n )) < e) > 1 - n. (26) 

n— too 

Now we note that the process (Z n ,Q n ) satisfies the conditions (a) compact containment prop- 
erty (Lemma 13. 2p and (b) oscillation bound ( ([25]) ) of Corollary 3.7.4 in [5] . Hence (Z n ,Q n ) is 
tight in D([0,oo),A4 x 1R + ) and (Z n ,Q n ) =^ (Z,Q) along some subsequence for some random 
variable (Z,Q) taking values in L>([0, oo), .A/f x K+). Also ([26]) satisfies the condition (a) in 
Theorem 3.10.2 in [5] which implies that (Z,Q) has continuous paths almost surely. 



4 Characterization of the limits 

In this section, we characterize some properties of the limits which lead to uniqueness. To have 
simple notations, we consider full sequence to converge instead of subsequence. To this end, 
we intend to define all the variable on a common probability space using Skorohod represen- 
tation theorem. From Condition 12.1( a). it is clear that E n is tight in D([0, oo), M + ). Hence 
(Z n ,Q n ,E n ) is tight in D([0,oo),M x M+) x D([0, oo),R+). Therefore by Skorohod repre- 
sentation theorem we can say that (Z ln ,Q ln ,E ln ,U ln ,V ln ) (2 1 ,Q 1 ,A-) in D([0,T],M x 
M+) x D([0, oo), R+) almost surely on some probability space (f2,.F,P) where 

law of (Z ln ,Q ln ,E ln ) = law of (Z n ,Q n ,E n ) for all n, 

and law of {Z , Q 1 ) = law of (Z, Q). 

Also (Z l ,Q l ) has continuous paths almost surely for Z\ = (1,-Zj). Hence for any T > 0, 
the followings hold, almost surely: 

= (27) 
= 0, 

= o, 

= 0, (28) 



lim 


sup 




z]\ 


n— >oo 


s£[0,T] 




lim 


sup 


\Ql n - 


Ql 


n— >oo 


s€[0,T] 






lim 


sup 


\z] n - 


z\ 


n— >oo 


SG[0,T] 






lim 


sup 


\B] n - 




n— >-oo 


se[o,T] 
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where Z\ = (1,2:/), B\ n = E\ n - Q\ n and B\ = Xt - Q\. Hence from © and ©, we get 
X ln -4 X 1 uniformly on [0,T],T > 0, and 

Xl = Q]+Z} and Q\ = [X} - 

(|28p implies that is nondecreasing in i and so it is a function of bounded variation on [0, T] 
for all T > 0. Also from above, it is easy to see that 

law of (Z 1 ,Q 1 ,5 1 ,X 1 )= law of (Z,Q,B,X). 

Lemma 4.5 Let (S, ir) be a metric space and K C S be compact. Let f : S — )• R be a function 
satisfying the following: for any sequence s n — > s and s £ if, /(s n ) — > /(s) as n-y oo. T/ien 
/or any e > 0, £/iere exists a 5 > snc/i i/toi |/(si) — /(s2)| < 6 whenever si £ S,S2 £ K and 
7r(si,s 2 ) < 5. 

Proof: If not, then there exists two sequences {s n }, {s„} such that {s n } C if and 

vr(s n , s n ) < -, |/(s n ) - f(s n )\ > e V n > 1. 
n 

Now if being compact, there exists s £ if such that along some subsequence {n^}, s nfc 4sG 
if. Hence s nfe — )• s £ If as rifc — >■ oo. This is contradicting to the fact that \f(s nk ) — f($n k )\ > e 
for all rifc. Hence the proof. □ 

Following lemma is a consequence of Lemma 13.31 

Lemma 4.6 Fix T > and xo £ R+. Consider a decreasing sequence {f n } in C&(R) so £/ia£ 
/ n > 0,/ n (x) = 1 on [xo — ^,#o + h) an d f n vanishes outside of [xq — ^,xq + jj]. T/ien 
IP(lim ri _ ! . 0O sup tg [ 0j T] (/", Z 1 ) > 0) = 0. In particular, for t > 0, 2/ /ias no atom ai xo almost 
surely. 

Proof: Let P^inin^oo sup tg [ ,T] (/ n > ^ 1 ) > K 2) > for some positive constant K2,T]. Then 
F(sup te[0iT] {f m ,Z}) > k 2 ) > r] for all m. Note that {Z\ : t £ [0,T]} is compact in M. Now 
from (j27|) and Lemma |4"31 we get, sup te [ 0) T] (/ m , ) — >■ sup tg [ 0j T] (/ m , 2 1 ) as n — > oo almost 
surely. Therefor using Fatou's lemma, for all m, 

liminf P n ( sup (f m ,Z?) > K2) = liminf P( sup {f m ,Z} n ) > rc 2 ) > P( sup {f m ,Z\) > k 2 ) > ??. 
te[o,T] te[o,T] te[o,T] 

Now we choose k > from Lemma 13.31 for e, 77 replaced by k 2 /2,?7/2. Thus if we choose m 
large enough, we get 

rj < liminf P n ( sup (f m ,Z?) > k 2 ) < liminf P n ( sup Z?[(x - 7) V0,x + j] > rc 2 ) < J, 
n ^°° *e[o,r] n ^°° te[o,T] 4 4 2 

which is a contradiction. This completes the proof. □. 

An immediate consequence of the above lemma is -2^([0, 00) = Zt((0, 00)) = Zt for all t > 0, 
almost surely. 



13 



Lemma 4.7 For any t>0,Z t satisfies the fluid model equation (k,X,u) given by 

Z t (C x ) = Z (C x + S(0,t)) + f G c (x + S(s,t))dB s , 

Jo 

almost surely where S(s,t) = J k(X u )du. 

Proof: It is enough to prove the above result for the process Z . For t > and C £ B([0, oo)), 
we have 

Z?(C) = Z (C + S n (0, t)) + I?(C), (29) 

where 

J-l , tj+i 



n 

j=0 i=B?.+l 



3 



for any partition {ijK = o of [0,t]. But Z\ n might not possess the same expression as Z™ as 
the stochastic variables v l n might not make sense on new probability space (Cl, J 7 , P). Define 
S ln (s,t) = J*k n (Xi n )du and S 1 (s,t) = J* h(Xl)du. Fix x G M.+ and let A T be a countable 
dense set in [0,T]. Then for any e > 0, 

P(sup \Zl n (C x )-2% n (C x + S ln (0,t)) - f G c (x + S ln (s,t))dBl n \ > e) 
*gA t jo 

= P n (sup \Z?(C x )-Z£(C x + S n (0,t))- f G c (x + S n (s,t))dBl l \>e) 
teA T Jo 

= P n ( sup \I?(C X ) - f G c (x + S n (s,t))dB^\ > e). (30) 
te[o,T] Jo 

Applying Lemma 13.21 for any positive 77, we have constant K such that P n (f2™ ) > 1 — 77 where 
S7y = {sup tg r T i X™ < K} for all n large. 

We choose 5 > and partitions {tj}]Zo such that maxi<j<.7(t) |ij — tj-i\ < 5 and 
su Pte[o,T] = < 00 • Recall that y n — >• as n — > 00 in probability. From Lemma f3.1l and 
Lemma 13.21 we choose M such that P n (fig) > 1 — 77 for all n large where f2g = {supr 0j yi B™ < 
f }. Hence on ng, (B t » +1 - A") < M for all n. Again for i G {Bg + + 2, . . ',B? j+ J, 

we have ij < r" < tj+i- Recall the event J7^(M, M) from (|B2p . For any e\ > 0, we have for 

< 3 < J(t) - 1, 

£ E £ G c (x + 5"(rr,t)) 

i=B? +1 

i 



4 = 5" +1 

< - A")(^(a + 5« fe+1 , < ))-G c (x + 5-(t J ,t))) + £l 
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on fi^(M, M) n fig, for all n large. Using the fact that p(v n ,v) — > as n — > oo , we get on 

n\(M, m) n fig n fi?, 



J-1 5 



E~ E S v?(C x+Sn{T n, t) )- [ G c (x + S n (s,t))dBl 



j=0 i=B"+l 
J-1 

^ E^+i " 5 S) + (^i+i.t) - ei) V 0) - G c (x + S n (^))) + Me 1 + J(t)e ls 

j=0 

< M\g\ 00 (K5 + e 1 ) + (M + J(5))e 1 , 

for all n large where \g\oo denote the supremum norm of g. First choosing 6 > small enough 
and then choosing e\ we can have the r.h.s. less than e/2 on Q^(M, MjOflgfin" for all n large 
and for all t G [0, T]. A similar calculation gives that ^"(C^) - Jjf G c {x + S n (s, i))ciB s n > -e/2 
on fi^(M, M) n fig n fi? for all n large and for all t G [0, T]. Since liminf^oo P„(fi^(M, M) n 
fis n fi") > 1 - 3r/, we have from (|50j) 

limsuplP( sup \Zj n {C x ) - Zl r \C x + 5 ln (0,t)) - f G c (x + 5 ln (s, t))d/5 s ln | > e) < 3r/. 

n->oo ie[0,T] JO 

7/ begin arbitrary, we have for any e > 0, 

limsu P P( sup \Z} n (C x ) - Zl n {C x + S ln {<d,t)) - I G c {x + S ln (s,t))dBl n \ > e) = 0. (31) 

n-s-oo te[0,T] JO 

Since sup [0>T] \X l s n - X l s \ -> 0, by Condition sup [0 T] |/c n (X s ln ) - &(Xj)| -> as n -> oo, 
almost surely. Hence 

sup \S ln (s,t) — S 1 (s,t)\ — > as n — > oo. 
se[o,T] 

This implies sup se[0 T] |G c (aH-S ln (s, t))-G c (x+S' 1 (s, t))| -»• as n -»• oo. Since /?(-2 ln , ^q 1 ) -»• 
almost surely, we have for t E [0, T] 

Zl n {C x +S ln {Q,t)) < Z^C x +S 1 (0,t)-e 2 )+e 2 ,Z^C x +S 1 (0,t)+e 2 ) < Z^ n (C x +S ln (0,t))+e 2 , 

for any chosen e 2 > and all n large (might depend on sample point). By Condition 12.1( c). 
Zq is deterministic with distribution function Lipschitz continuous and so sup tg [ 0] T] \^o n {C x + 
S ln (0, t))-_Zl(C x +S l {Q, t))\ ->■ almost surely. From PSJ), it is easy to check that p{dB ln , dB 1 ) 
where dB ln ,dB l are considered as Borel measures on [0,T]. Since B 1 is continuous almost 
surely, applying Theorem A2.3.I in [3] and (|28p . we get 

sup | / G c (x + S 1 (s,t))dBl n - [ G c (x + S 1 (s,t))dBl\ ^ as n oo, almost surely, 
te[o,T] Jo io 

and hence 

sup | f G c {x + S ln (s,t))dBl n - [ G c (x + S 1 {s,t))dB 1 s \ ^0 as n ^ oo, almost surely. 
te[o,T] jo jo 
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Now we show that sup 46 [ 0j T] \Z^ n (C x ) — Zj(C x )\ — > as n — > oo almost surely. Consider the 
map / : Ai — > R defined by f(fJ.) = ^(C x ). From Lemma 14.61 we see that / satisfies the 
condition of Lemma 14.51 for the compact set {Z^ : t G [0,T]}, almost surely. Hence using 
([271) and Lemma we have sup fg [ 0T ] \Z\ n (Cx) — Zl(C x )\ — > as n — > oo, almost surely. 
Combining the above estimates with (|3ip . we get Q(x,T) G T such that P(f2(x,T)) = 1 and 

Z\(C X ) = Z^(C x + S 1 (0,t)) + f G^x + S^tydBl, 

Jo 

for all t G [0, T] on S7(x,T). Since {C x : x G R + ,:e rational} determines any Borel-measure 
uniquely on ffi_|_ we can take r^oo — n^gpj ^{x-.x rational} Q(x, T) on which 

Z\(C X ) = Z^(C x + S 1 (0,t))+ I G^x + S^s^dBl 

Jo 

for all t > and x G R + . This completes the proof as P(0 oo ) = 1. □ 

5 Appendix 
5.1 Appendix A 

In this section, we prove existence of unique solution to the fluid model type equations. 

Lemma 5.8 Fix T > 0. Let k, Hi : R — > R, i = 1,2,3,4, be bounded, Lipschitz continuous. 
Then the following integral equation 

ft ft ft ft ft 

x t = #i(/ k(x s )ds) + H 2 ( k(x u )du)ds + H 3 (x s )H 4 ( k(x u )du)ds (32) 

Jo Jo Js Jo Js 

XQ = ffl(0). 

has a unique solution in C([0, T],R). 

Proof: To simplify the notation, we define S(s,t,4>) = J"* k((f)(u))du for (j) G C([s,i],R). As- 
sume that solution x t is uniquely defined on [0,io] for to G [0,7"). We consider the following 
integral equation for t G [to,T] 



fto 

x t = Hx(S(Q,t ,x) + S(t ,t,x)) + / H 2 (S(s,t Q ,x) + S(t ,t,x))ds 

Jo 

rto rt 

+ / H 3 (x s )H 4 (S(s,t ,x) + S{t ,t,x))ds + H 2 {S{s,t,x))ds 

Jo Jt 

+ f H 3 {x s )H 4 (S(s,t,x))ds. (33) 

Jt 

Now define a operator F : C([t ,T],R) C([t ,T],R) as follows: 

= H 1 (S(0,t ,x) + S(t ,t,4>))+ H 2 (S(s,t ,x) + S(t ,t,^))ds 

Jo 

+ [ ° H 3 (x s )H 4 (S(s,t ,x) + S(t ,t,<j)))ds + [ H 2 (S(s,t,<j)))ds 
Jo J to 



+ / H 3 (cf>(s))H 4 (S(s,t,<j>))ds. 



to 
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To simplify the notation, we denote the i-th term on the r.h.s. of the above expression by 
Fi((fi) for i = 1,2,3,4,5. We denote the supremum (Lipschitz constant) of Hi by Hi^Li) for 
i = 1,2,3,4. Let L k be the Lipschitz constant of k(-). Then for <f> ,<jr G C([tQ,T],M) the 
followings hold: for t G [to,T] 





-Fi{ 


4> 2 )(t)\ 


< 


\m l ){t) 


-F 2 ( 


f)(t)\ 


< 




-F 3 ( 


<t> 2 )(t)\ 


< 




-m 2 )(t)\ 


< 




-F 5 ( 


4> 2 ){t)\ 


< 



hooLiLkt^t - to)\4> 1 — (j) 2 \t t, 
2 L k {t - 1 ) 2 \4> 1 -<p 2 \t t, 

hocUL k (t - i ) V - 4>% t + H ioo L 3 {t - to)^ 1 - 4> 2 \t t- 

Hence combining the above expressions we get, for t G [icb^l 

\F{(t> l ){t) - F(ct> 2 )(t)\ < (L\L k + L 2 L k T + H 3oo L A L k T + H Aoo L 3 )(t - t )\<f, x - <j) 2 \ tot . 

Hence we can choose h > small enough so that supj to t ] \F(c/) 1 )(s) — F((j) 2 )(s)\ < gift 1 — <j) 2 \t t 
for some positive g < 1 and t — to = h. So by contraction mapping theorem, there exists a 
unique continuous function x defined on [to , t] satisfying (i33j) . 

Putting to — 0, we see that xt satisfies ([32]) on [Q,h]. Having the solution defined on 
[0, nh AT], we can extend it uniquely on [0, (n + l)h A T] for n G N. Since h > is fixed, this 
defines the solution uniquely on [0, T]. □. 

We can extend Lemma 15.81 as follows: 

Lemma 5.9 Fix T > 0. Let k, Hi : R —> R, i = 1,2,3,4,5, be Lipschitz continuous. We also 
assume H A ,H§ to be bounded. Then the following integral equation 

x t = H 1 {S{0,t)) + [ H 2 {S{s,t))ds + [ H 3 {x s )H 5 (x s )H 4 (S(s,t))ds (34) 
Jo Jo 

xo = #i(0). 

has a unique solution in C([0, T],R) where S(s,t) = k(x u )du. 

Proof: Let ip n : R — > R be a smooth cut-off function such that < (p n < 1, (f n (s) = 1 on [— n, n] 
and (f n (s) = outside of [— n— 1, n + 1]. Define H™(s) = ip n (s)Hi(s) for i = 1, 2, 3. Then Hf is 
a bounded, Lipschitz continuous function for i = 1,2,3. Hence H^H§ is a bounded, Lipschitz 
continuous function for n G N. Therefore applying Lemma f5.8| we have unique x n : [0, T] — > R, 
continuous, satisfying 



xV 



f k{x n s )ds) + / flj( /"* k(x™)du)ds 

Jo Jo Js 

+ f H$W)H 5 (x2)H 4 ( f k{x n u )du)ds (35) 

J0 



x o = fl"i(0). 
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Since k,Hi,i = 1, ... ,5 are Lipschitz and H^,H§ are bounded, we can get positive constants 
d\ , d 2 such that 

\Hiif H^)ds)\ < di + d 2 ( \x n s \ds 
Jo Jo 

\[ H%(x n s )Hs{x n s )H±{ J k{xl)du)ds\ < d! + d 2 [ \x^\ds, 

Jo Js Jo 

and 

1/ H%( k{xl)du)ds\ < di + d 2 / / \x 7 l\duds<d l +d 2 T j \x n s \ds, 

Jo Js Jo Js Jo 

for all t G [0, T] and n G N. Now combining these estimates with (135p and applying Gronwall's 
inequality we have 

sup sup | x™ | < d$ , (36) 

n [o,T] 

for some constant d%. For any compact CcK, there exists constant Lq such that \H™(x) — 
H'i{y)\ < L c \x — y\ for all x,y G C and n G N. Therefore using an analogous expression as 
(|33|) . we get a constant c?4 > (depending on d%) satisfying 

sup sup |x™ — x™| < di\t — s\. (37) 

n 0<s<t<T 

(I36D and (I37p imply that the sequence {x n } is equi-continuous family of continuous functions 
on [0, T]. Therefore using Arzela-Ascoli, theorem there is a x : [0,T] — > K, continuous, such 
that \x nk — x\ot — > along some subsequence nk — > 00. Hence letting — > 00 in (j35j) we have 

x t = Hi( J k(x s )ds) + / H 2 ( I k(x u )du)ds 
Jo Jo Js 

+ / H s (x s )H 5 (x s )H A ( / k(x u )du)ds 
Jo Js 
xo = x = H 1 (0). 

This proves the existence of solution to ()34p . To prove uniqueness, let x be another solution 
to (I34p . Define a n = inf{t > : \xt\ > n} A T. Since H n (s) = H(s) for \s\ < n, from Lemma 
15.81 we get x s = x™ = x s for s < cr n and for all n large (we need to take large n to ensure that 
xo G [— n,n]). Therefore to complete the proof it is enough to show that lim inf n _s.oo &n — T. 
But this is obvious as sup[ 0i T] \%s\ < 00 (follows from a simple calculation similar to ([36]) ). □ 

5.2 Appendix B 

Consider a sequence of probability measures {v 11 } and v on [0, 00) such that v n — > v as n — > 00. 
Let {vf J-^ffoo be an i.i.d. sequence with common probability distribution v n . For m G Z and 
4 > 0, define 

m+|n€J 

£™(m,£) = - <V ( B1 ) 

i=l+m 
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By Skorohod representation theorem, there exists [0, oo)-valued random variables Y n ~ v n 
and Y ~ v such that Y n — > Y almost surely on some common probability space. Define 
Y = sup n Y n . Let v be the law of Y. There exists a continuous, increasing, unbounded 
function / such that / > 1 and (f 2 ,v) < oo (see Appendix B in |18j). Define 

V = { XCx ,x > 0} U {xc x ,x > 0} U {/"}. 

Lemma 5.10 Fix M,L > 0. Under the above assumptions, for all e, r/ > we have 

limsupP n ( max sup sup | (/, C n (m, £)) - £(f, v n ) \ > e) < 7]. 

n-^oo - nM<m<nM ££[0,L] feV 

For the proof of above lemma we refer Lemma B.l in |18j . Following the same argument as 
in Lemma 5.1 in [TB] we can have a sequence e^n) such that — > as n — > oo and 
lim^oo P n (r2^(M, L)) = 1 for every fixed M, L > where 

Q n A (M,L) = { max sup sup\(f,£ n (m,£))-£(f,v n )\ <e A (n)}. (B2) 

-nM<m<nM te [ 0>L ] / eV 

Acknowledgement: The author is grateful to Prof. Rami Atar for his valuable suggestions 
to improve this paper. 
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